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Abstract. We analysed the optical and infrared properties of X-ray sources in the 
Great Observatories Origins Deep Survey (GOODS), a deep, multiwavelength survey 
covering 0.1 square degrees in two fields. The HST ACS data 1 are well explained by a 
unified AGN scheme that postulates roughly 3 times as many obscured as unobscured 
AGN, as are the spectroscopic and photometric redshift distributions once selection 
effects are considered. Our model predicts infrared number counts of AGN that agree 
well with the preliminary Spitzer data, confirming that large numbers of obscured AGN 
are present in the early Universe (z > 1). 



1 Why There Might be Obscured AGN at High Redshift 

Luminous Active Galactic Nuclei (AGN) have been readily found throughout 
the Universe, primarily in surveys at radio, optical, ultraviolet, and X-ray wave- 
lengths. These AGN (often called quasars) are ultimately powered by gravity, 
as matter from the host galaxy accretes onto a central supermassive black hole. 
The prevalence of AGN at redshifts z ~ 2 — 3 indicates that that is the epoch of 
greatest black hole growth. It is also an epoch of strong star formation in young 
galaxies, one of several arguments that the growth of supermassive black holes 
must be closely tied to the formation of galaxies. Understanding the co-evolution 
of galaxies and black holes requires knowing the accretion history of supermas- 
sive black holes, via an accurate census of black hole demographics: how many 
there are, where they are, and what their masses are. 

Although many thousands of AGN have been catalogued, most recently in 
the Sloan Digital Sky Survey (e.g., Richards et al., 2002; Anderson et al., 2003), 
the demographics of supermassive black holes has proved a surprisingly elusive 
goal. Very few "type 2 quasars" 2 have been found (e.g., Norman et al., 2002; 
Stern et al., 2002; Dawson et al., 2003). Of course, UV-excess or optical emission- 
line surveys would not have found most obscured AGN, nor would soft X-ray 
surveys such as the ROSAT All-Sky (Voges et al., 1999) or WGA (Singh et al, 

1 Based on observations obtained with the NASA/ESA Hubble Space Telescope, which 
is operated by the Association of Universities for Research in Astronomy (AURA), 
Inc., under NASA contract NAS5-26555. 

2 "Type 2 quasars" are the luminous analogs of local Seyfert 2 galaxies, i.e., AGN 
lacking broad optical emission lines. According to the unification hypothesis for AGN, 
the central continuum and broad line region in type 2 AGN are obscured by high 
column densities of gas and dust (Nh > 10 22 cm -2 ) along the line of sight. 
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1995) surveys. Instead, one needs to look at hard X-rays, where absorption is a 
smaller effect, or in the far-infrared, where the absorbed energy is re-radiated. 

There are three good reasons to expect many supermassive black holes to 
be hidden behind a thick screen of gas and dust. First, local AGN conform 
well to the unification scenario, in which roughly 3/4 of all AGN are at least 
partly hidden from a direct line-of-sight (Antonucci, 1993). Second, the 2-40 keV 
spectrum of the X-ray background is very hard, much harder than the typical 
spectrum of an unobscured AGN (Mushotzky et al., 1993), yet high resolution X- 
ray imaging makes clear that the "background" is in fact summed emission from 
individual AGN (Giacconi et al., 1979). As X-ray astronomers have understood 
for nearly two decades, this means most AGN are absorbed by column densities 
in excess of 10 22 atoms/cm 2 (Setti & Woltjer, 1989; Madau, Ghisellini, & Fabian, 
1994; Comastri ct al., 1995; Gilli et al., 2001). Recent papers (Comastri et al., 
2003; Worsley et al., 2004; Ueda et al., 2003; U03; see also Fabian, this volume) 
have argued that to explain the ~ 30 keV peak in the energy density of the X-ray 
background requires a substantial contribution from Compton-thick AGN (with 
obscuring column densities > 10 24 cm -2 ). Third, in the young Universe, where 
forming galaxies are often dusty and reddened, obscuration of the central AGN 
is even more likely. Thus, at traditional survey wavelengths, the strong line and 
nonthermal continuum emission by which AGN are selected and identified may 
be completely invisible. 

Our strong reliance on optical identification has led many to argue that 
luminous obscured AGN, few of which have been found at high redshift, do 
not exist, and the obscuring material must be blown away in higher luminosity 
sources. While this possibility cannot be ruled out, it is certainly premature 
to explain the absence of obscured AGN from surveys that were remarkably 
insensitive to such objects. 

2 Searching for Obscured AGN in GOODS 

A strong motivation for the Great Observatories Origins Deep Survey (GOODS) 
was the discovery of obscured AGN. GOODS consists of deep imaging in the far 
infrared with the Spitzcr Space Telescope (Dickinson & Giavalisco, 2002) and 
in the optical with the Hubble Space Telescope (Giavalisco et al. , 2004) on 
the footprints of the two deepest Chandra fields (Giacconi et al., 2001; Brandt 
et al., 2001; Alexander ct al. , 2003). The total area is roughly 60 times larger 
than the original Hubble Deep Field (Williams et al., 1996) and nearly as deep in 
the optical. The Great Observatories data were augmented with ground-based 
imaging and spectroscopy. 3 More than 85% of the GOODS X-ray sources are 
AGN, with luminosities Lx > 10 42 ergs/s; to avoid confusion with starburst 
galaxies, we limit our sample to X-ray sources above this luminosity threshold. 

With extensive coverage over five decades in energy from 24 |im to 8 keV, the 
GOODS survey is well suited to finding a high-rcdshift population of obscured 

3 GOODS observations are summarized at http://www.stsci.edu/science/goods/. 
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Fig. 1. X-ray fluxes and z-band magnitudes of GOODS AGN. (Left) Hard X-ray flux 
distribution for X-ray sources in the combined GOODS-North and -South fields (heavy 
solid line), compared to number counts calculated from a simple unification model 
(light solid line). Contributions from unobscurcd (dashed line) and obscured (dotted 
line) AGN are shown separately; these follow similar distributions because hard X- 
rays are not strongly affected by absorption. (Right) Distribution of observed z-band 
magnitudes for GOODS-North and -South X-ray sources (heavy solid line), compared 
to model distribution (light solid line). Obscured AGN (dotted line) dominate at faint 
magnitudes, where only the host galaxy contributes to the optical light, while bright 
samples consist both unobscured (dashed line) and obscured AGN. 



AGN. A complementary approach, given the relatively low surface density of 
AGN compared to normal galaxies, is to target higher luminosity AGN over 
a wider area of the sky, an approach followed by the ChaMP (Green et al., 
2004), CYDER (Castander et al., 2003), and HELLAS2XMM (Fiore et al., 2003) 
surveys, among others. Here we describe the AGN detected in the X-ray and 
optical in the GOODS North and South fields, which represent an order of 
magnitude more objects than in most previous works. 

Most studies of black hole demographics begin with the AGN content of a 
given survey, correcting where possible for selection biases to infer the underlying 
population. If selection effects are strong, however, one ends up making large 
extrapolations using little information. We therefore took a different approach: 
we asked, if there is a substantial population of obscured AGN, what would be 
seen in a deep multiwavelength survey like GOODS? This work is presented by 
Treister et al. (2004a); here we briefly describe our assumptions and results. 

We start with a very simple unification scheme for AGN, in which a cen- 
tral engine - well understood from studies of unobscured AGN - is surrounded 
by a dusty torus 4 that absorbs optical through soft X radiation and re-radiates 

4 Note that the geometry need not be a torus; we use this geometry simply to derive 
the distribution of absorbing column densities expected from random orientations. 
The resulting Nh distribution matches well the observed one, with the proviso that 
Compton-thick sources are missing, but so would many other geometries. 
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the absorbed energy as thermal infrared emission. We combine well-measured 
spectra of unobscured AGN with plausible reddening and absorption laws, and 
with the dust emission models of Elitzur et al. (2003) and Nenkova et al. (2002), 
to create a grid of spectra appropriate to increasing absorbing column densities 
along the line of sight to the AGN. The torus geometry was adjusted so that 
three-quarters of all AGN are obscured (i.e., have Nh > 10 22 cm -2 ) and so 
the mid-plane column density is 10 24 cm~ 2 (even the Chandra deep surveys are 
not sensitive to objects with higher column densities). We further assumed this 
geometry was independent of redshift and luminosity - the simplest assumption, 
and one that, frankly, we expected to prove wrong. We added an L* host galaxy 
spectrum to each AGN template spectrum. Finally, we assumed the most re- 
cent, hard- X-ray-selected luminosity function and luminosity-dependent density 
evolution given by Ueda et al. (2003; U03). 

As shown in Figure 1, this very simple model explains very well the X-ray 
and optical counts in the GOODS fields. That it explains the X-ray counts is not 
too surprising since we start with a hard X-ray luminosity function and evolution 
derived from these fields, among others. For the optical, however, the agreement 
indicates that faint optical sources are primarily obscured AGN. Indeed, these 
tend to have red colors and hard X-ray spectra, as expected of obscured AGN. 

These results mean that the GOODS HST data agree well with a popula- 
tion of obscured AGN extending from z = to well beyond z > 1. Specifically, 
the X-ray and optical properties of GOODS AGN are fully consistent with a 
very simple unification picture, requiring no modification with cither redshift or 
AGN luminosity. The implication is that some X-ray-emitting AGN are missed 
in the Chandra Deep Field samples, and others lack optical spectroscopic iden- 
tifications. Since the best 8- or 10-m class telescopes have an effective limit at 
R < 24 mag, many obscured AGN will lack spectroscopic identification. 

The redshift distribution is a critical point. The best population synthesis 
models for the X-ray background prior to the new Chandra and XMM deep 
fields (Gilli et al., 2001) incorporated a high-redshift population of obscured AGN 
that evolves similarly to the unobscured AGN, and hence predicted a redshift 
distribution peaking at z ~ 1.5. Yet observed redshift distributions from the 
Chandra Deep Fields and the Lockman hole peak at much lower rcdshifts (e.g., 
Hasinger, 2002; Barger et al., 2003). This has been taken as evidence of a fatal 
flaw in the population synthesis models. 

We find the discrepancy does not imply any lack of obscured quasars at high 
redshift. Rather, it results from a combination of two factors: (1) the luminosity- 
dependent density evolution now measured for X-ray AGN, which was not used in 
the Gilli et al. (2001) models and which increases the fraction of AGN at lower lu- 
minosity and lower redshift, and more importantly, (2) the implicit spectroscopic 
limit for X-ray samples, which preferentially excludes the faint counterparts of 
high-redshift, obscured AGN. 

The latter point is illustrated in Figure 2, which shows the redshift distri- 
bution for GOODS AGN. The upper dashed line shows the predicted redshift 
distribution for our unified model, which still peaks above a redshift of 1, much 
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Fig. 2. Observed redshift distributions (heavy solid line) for AGN in the GOODS-South 
(left panel) and -North (right panel) fields, including both spectroscopic (hatched area; 
Barger et al., 2003; Szokoly et al., 2004) and photometric redshifts (Barger et al., 
2003; Mobasher et al., 2004); these are 100% complete in the GOODS-S field and 75% 
complete in the GOODS-N region. The predicted distribution {upper dashed line) peaks 
above a redshift of 1, similar to earlier population synthesis models and in apparent 
disagreement with the observed one. However, if we cut the sample at R < 24 mag 
(lower dashed line), as appropriate for a sample with optical spectra, the model and 
observed distributions agree well. The data are therefore consistent with a significant 
high-redshift population of obscured AGN which are missed in spectroscopic samples 
due to their faint optical magnitudes. 



like the earlier population synthesis models. If we impose a spectroscopic limit 
of R = 24 mag, however, we obtain the lower dashed line, which agrees well 
with the observed spectroscopic redshift distribution. Note that the photometric 
redshifts are indeed shifted to higher redshifts compared to the spectroscopic 
redshifts, although they too miss the very faintest, highest redshift obscured 
AGN. 



3 Spitzer Data and the Far-Infrared Properties of AGN 

According to the unification scheme, dust and gas absorbs much of the luminous 
optical through X-ray AGN emission, and is thus heated. The absorbed energy 
is re-radiated as thermal emission, which should be readily visible in the far- 
infrared. A key element of the GOODS project is very deep imaging with the 
Spitzer observatory, in the IRAC bands from 3.6 to 8.4 microns and the MIPS 
24-micron band. These data will provide a strong test of our unified model. 

We used our simple unification model to predict the far-infrared source 
counts, in each Spitzer band, for the appropriate GOODS flux limits (Treis- 
ter et al., 2004a). As of spring 2004, some IRAC data in the GOODS-S field 
were already in hand. We compared preliminary photometry of the AGN to our 
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Fig. 3. Predicted Spitzer counts for the simple unified model. (Left) The model pre- 
diction at 3.6-u.m (dashed line) agrees well with the preliminary Spitzer counts of 
GOODS-S X-ray sources (solid histogram), in both normalization and shape. (Right) 
Model prediction at 24 urn for total source counts, (upper dot-dash line) and for AGN 
brighter than the X-ray flux limit of GOODS (solid line); the difference represents 
the AGN missed by Chandra. Dotted and dashed lines represent the obscured and 
unobscured AGN, respectively; the former dominate by the specified factor 3:1. 

source count prediction. The results are shown in Figure 3 (van Duyne et al., 
in prep.), along with predictions for the 24-micron source counts. (The first 
GOODS MIPS data will be obtained in summer/fall 2004.) 

Certainly these early data require closer inspection, and the model awaits 
a stronger test from the full multiwavelength Spitzer data set. The pending 
data will allow us to verify whether large numbers of obscured AGN are indeed 
present in the early Universe, and will also be extremely important for refining 
the AGN model, specifically the covering factor, dust geometry, and evolution 
thereof. The 24-micron images are particularly powerful, since the AGN emission 
is most isotropic at the longest infrared wavelengths. Obscured AGN should be 
very bright far-infrared sources; those that are missed even by X-ray observations 
probably look like ultraluminous infrared galaxies. 

We have looked at the spectral energy distributions of the GOODS AGN, 
including the available Spitzer data. Indeed, the hard X-ray sources with faint, 
red optical counterparts are extra-luminous in the Spitzer bands. Two example 
spectral energy distributions are shown in Figure 4 (van Duyne et al., in prep.). 
Both cases look just like a normal galaxies at optical wavelengths, but are re- 
vealed to be AGN by the luminous X-ray emission. They are clearly obscured by 
virtue of their faint optical magnitudes, hard X-ray spectra, and bright infrared 
fluxes. 
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Fig. 4. Spectral energy distributions of two GOODS AGN, including HST ACS data 
(squares), Chandra X-ray data (triangles, straight lines), and Spitzer data (filled cir- 
cles). Also shown are the best-fit Seyfert 2 galaxy templates, both observed (light solid 
line) and model spectra (heavy solid line). 



4 Summary and Conclusions 

The GOODS multiwavelength survey supports the idea of a significant popula- 
tion of obscured AGN persisting to high redshift. From a strictly observational 
viewpoint, we cannot rule out this population, and indeed, sec every evidence 
for its existence. In a subsequent paper (Treistcr ct al., 2004b) we explain the 
apparent decline in the ratio of obscured to unobscured AGN with redshift or 
luminosity, which is seen in several X-ray-selected AGN samples (Hasinger, 2004; 
U03), by the same selection effects against Compton thick and/or optically faint 
sources. 

The data are consistent with a simple unified model of a luminous AGN nu- 
cleus obscured along some lines of sight by dust and gas. Our work is not strongly 
sensitive to the geometry, and other configurations of dust are equally possible as 
long as they predict a similar distribution of absorbing column densities. We are 
also not sensitive to the numbers of Compton-thick AGN, a population that will 
be probed with future X-ray missions sensitive in the energy range 10-100 keV. 
In any case, we conclude that at least 3/4 of all supermassive black holes are 
obscured by column densities of at least Nh ~ 10 22 cm~ 2 . 

The unified model implies that half the obscured AGN are not detected in 
X-ray fields at the Chandra Deep Field depths. They should be clearly visible in 
the Spitzer data, although separating heavily obscured AGN from very luminous 
starbursts will not be simple. Upcoming Spitzer data, along with the integral 
constraint from the X-ray "background" spectrum, will better constrain the 
geometry and other properties of the obscuring material. 

Individual spectral energy distributions of the AGN will also constrain the 
dust geometry, and perhaps the contribution from any starburst component. 
Bolometric corrections in the far- infrared are likely to be significant, and affect 
estimates of the mass accretion rate and thus the growth rate of supermassive 
black holes. Note that observed spectral energy distributions never represent 
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the actual output from the AGN engine: the luminosity of unobscured AGN is 
overestimated (since the optical through soft X-ray emission is not isotropic), 
while obscured AGN appear under-luminous relative to their intrinsic power. 
Bolometric corrections can be calculated from models, such as those considered 
here, but can not be derived empirically for any highly anisotropic source. This 
is critical for understanding the energetics and structures of AGN. 

One strong conclusion from our work is that one should look for obscured 
black holes where they can be easily found, namely in the far-infrared and hard 
X-ray. Any substantial and widespread obscuration has important implications 
for the cosmic accretion history of black holes, and for the average efficiency of 
converting gravitational potential energy to radiation. 
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